An Effective Field Model to High Temperature Superconductors: 
Fitting of The Energy Gap to the Experimental Results of the 

Bi^SriC 'aCu20s 

E.C. BastoneQand P.R. Silva 

Departamento de Fisica, Instituto de Ciencias Exatas, 

Umversidade Federal de Minas Gerais, 
Belo Horizonte, CEP 30.123-970, CP. 702, MG, Brazil 
e-mail: erika@fisica.ufmg.br 



ABSTRACT 

We use an effective field model (transverse fsing model) to describe the de- 
pendence on the temperature of the energy gap of some two-dimensional (2 — D) 
superconducting systems. The order parameter of this model is put in a direct cor- 
respondence with that of the Ising model. Then, we use the exact relations for the 
spontaneous magnetization of the fsing model in some 2 — D lattices as a means to 
fit the experimental results of Bi^Sr^C aCu^O^ films obtained through reflectivity 
measurements by Brunei et al [Phys. Rev. Lett. 66, f346 (f99f)] and the argon- 
annealed Bi2Sr2CaCu20s samples investigated by Staufer et al [Phys. Rev. Lett. 
68, f069 (f992)] through Raman scattering. The zero temperature energy gap E g (0) 
relations were evaluated in various 2 — D lattices and also in the (3 — D) simple cu- 
bic lattice, for comparison. In the square lattice case we obtained E g (0) — 3.52fcT c , 
coincidentally the same number as the BCS result, and in the triangular lattice case 
we get E g (0) = 3.30fcT c in agreement with the experimental findings of Brunei et 
al. 



Pacs number: 74.20.De, 74.72.Hs 

As has been pointed out by Anderson 1 : the consensus is that there is absolutely no consensus 
on the theory of high temperature superconductivity. Anderson 2 attributes the novel phenomenol- 
ogy present on cuprates materials to a second kind of metallic state, namely, the Luttinger liquid. 
Schrieffer 1 has pursued the interplay between antiferromagnetism and superconductivity, extend- 
ing the BCS pairing theory beyond the Fermi- liquid regime in terms of spin polarons or " bags" . 
As pointed out by Cox and Maple 3 , the cuprates materials possess strong electronic correlations, 
where by " strong correlations" we mean that the average interaction energy substantially exceeds 
the average kinetic energy of the partially filled 3c? state of the Cu 2+ ion. Therefore the discussion 
of what would be the proper mechanism to describe the high temperature superconductivity in 
cuprates becomes polarized between the BCS theory defenders as Schrieffer, and the people who 
trust in the Hubbard model description as Anderson 1 .This controversy suggest us that the propo- 
sition of new effective hamiltonians, which does not enter in specific details of the mechanism, could 
be an alternative way to look at this problem. Inspired in the Gorter-Casimir two-fluid model of 
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superconductivity 4 , and also taking into account the Blinc-de Gennes model 5 ' 6 for the hydrogen- 
bonded ferroelectrics, Gaona J. and Silva 7 have proposed the following effective Hamiltonian as a 
means to describe some basic features of the superconducting state: 

H=-il^X i -^J2 J vZ i Z j . (1) 

The second term of (1) is a Ising-like term, where the operator Z is related to the wave function of 
the condensate of holes (electrons) and the coupling Jy favors the pairing of holes at two different 
lattice sites, contributing to the coherence of their wave functions below the critical temperature 
T c . In the first term of (1) the operator X is related to the wave function of the normal electrons 
(holes), where Q represents the transverse field. This term, which gives dynamics to the model, 
could represent a mimic for the motion of the free electrons (holes) through the Fermi barrier, 
in close analogy with the tunneling of protons through the double-well potential barrier in the 
KH2P04-likc (hydrogen-bonded ferroelectrics) case (see 5 ). 

We treat the dynamics of Hamiltonian (1), by applying the Random Phase Approximation 
(RPA), where we replace a time dependent expectation value < Si > t by a constant part < Si > 
- which is just the MFA expectation value - plus a small time dependent deviation S < Si > e' iut 
from the molecular field solution, to the Heisenberg equation of motion. The solution of it in the 
first order approximation is 7 : 

lo 2 = J 2 <Z > 2 +n (fi - J < A >) . (2) 

This solution is constrained by the zero order approximation 

(n - j < x >) < z >= o . (3) 

Below T c , < Z >^ and by putting (3) into(2), leads to 

u)± = ± J < Z > . (4) 

In (4) lo± can represent a two level system, and we could treat it as a gas of elementary excitations 
related to the hole (electron) condensate. It seems natural to interpret the separation in energy of 
this effective two- level system as the energy gap E g (T). We have (in units where h = 1) 

2A (T) = E g (T) = uu- (T) - u+ (T) = 2 J < Z > . (5) 

Now, by using an idea introduced by Rummer and Ryvkin 8 , since the quasi-particles are non 
interacting fcrmions and their numbers are not fixed, then the averaged occupancy of the levels 
u>± (T) is determined by the Fermi-Dirac distribution formula with zero chemical potential. Then 
we have 

<n(w)>= [e^ + 1]" 1 , (6) 

where [3 — 1/kT. By using the above distribution and associating a |— } eigenstate to the uj+ energy 
and a |+) eigenstate to the w_ energy, we obtain the following relation for the superconductor order 
parameter < Z >, namely 

< Z >= (-1) (n+> + (+1) (n_) , (7) 

or alternatively we can write 

< Z >= tanh pJ < Z >) /2] . (8) 
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Close to T c , < Z > goes to zero and we have from (8): 



2kT c = J . (9) 

At T = 0, < Z >= 1 and by putting (9) into (5) we obtain the energy gap 

E g (0) = 2A (0) = 2 J = AkT c . (10) 

It is interesting to mention that Hao 9 have obtained this value for the energy gap as one 
of the possible solutions of the BCS 10 self-consistency equation. Equation (8) can be put in 
a direct correspondence with the self-consistency equation for the Ising model in the mean field 
approximation (MFA), but with the temperature multiplied by a factor of 2 11 . Now let us consider 
the situation where the thermal correlations close to T c are very important, so that the MFA is a 
bad description to the order parameter behavior of the system. Indeed this seems to be the case 
in some two-dimensional superconducting systems, as for instance the Bi 2 Sr 2 CaCu20s, where 
the order parameter were experimentally determined through Raman scattering and reflectivity 
measurements 12 ' 13 . The temperature dependence of the order parameter of the Ising model 
in two-dimensions is exactly known 14 . We can assume that, for a two-dimensional (2 — D) 
superconducting system where thermal correlations can not be neglected, the order parameter 
can be mapped into the exact relation for the spontaneous magnetization of the 2 — D Ising 
ferromagnetic model. However in this case and in the neighborhood of T c , the gas of elementary 
excitations approximation does not work quite well and we must turn to the Boltzmann statistics. 
Now let us consider the exact relation for the spontaneous magnetization of the Ising model in the 
square lattice 14 ~ 16 : 



m 



2tanh 2 (2/?J)-l l] 



tanh 4 {213 J) 



(11) 



In figures 1, 2 and 3, we compare eq. (11) with the experimentally determined order parameter 
in BiiSriCaCuiO% single crystals, investigated through Raman scattering by Staufer et al. 12 
and in films of the same material obtained through reflectivity measurements by Brunei et al. 13 . 
For comparison, we plot in the same figure the eq. (8) and the BCS curve 10 . We observe that 
the curve given by eq. (11) nicely fits to the measurements of the energy gap performed at the 
Bi 2 Sr 2 CaCu20$ films, a manifestly 2 — D system. On the other hand, in the case of single crystals, 
the adjust to curve (11) is good in the case of the essentially 2 — D argon-annealed samples. In 
the oxygen annealed crystals the energy gap can be described both by BCS theory as equivalently 
by eq. (8). In fitting the experimental points to eq. (11), we have considered that: 

m -^(0)-2A(0)' (12) 

and in the case of eq. (8) m coincides with < Z >. 

Now let us consider the exact result for the critical temperature of the Ising model in the square 
lattice 14 < 15 . We have 

^ = 2.27 (13) 
which leads to the zero temperature energy gap 

HkT 

2A (T) = 2 Jo = 2zJ = = 3.52kT c , (14) 
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where we took z = 4 (the coordination number of the square lattice). The result we got in (14), is 
coincidentally the same as the BCS result and must be compared with the experimental finding of 
3.3kT c , obtained by Brunei et al. 13 . It must be emphasized that although the zero temperature 
energy gap has been derived from a mean field analysis (eq. (5)), its use in conjunction with the 
exact result for T c (eq.(13)) does not lead to contradictions, once at T = 0, we are far from the 
critical temperature T c , where thermal correlations are important and the mean field treatment 
substantially differs from the exact one. 

Until now, we have compared the experimental measurements of the order parameter of the 
Bi2Sr2CaCu20s with the exact result of the spontaneous magnetization of the Ising ferromagnetic 
in the square lattice. It would be interesting to pursue further on this subject, by comparing within 
the same token, the zero temperature energy gap obtained by considering other 2 — D and 3-D 
lattices. Therefore in the following we present the relations for the spontaneous magnetization 
(SM) of the Ising model in some other 2 — D lattices. 

In the case of the honeycomb lattice 17,18 we have: 



to 



1 - 2y / l-tanh 2 (2/3 J) 



tanh 2 (2(3 J) 



2 - tanh 2 (2/3 J) - 2 tanh 2 (2/3 J) - 2yT - tanh 2 (2/3 J) 



1/8 



(15) 



For the triangular lattice case, the order parameter behavior is given by the equation 



15,19. 



TO : 



2 tanh (2/3 J) - 1 



2 tanh 3 (2/3 J) - tanh 4 (2/3 J) 



1/8 



(16) 



In the case of the simple cubic lattice we do not have an exact relation for the SM. However, fifteen 
years ago, one of the present authors 20 was able to "guess" a relation for the SM of this 3-D 
lattice, which reproduces very closely the series results of Fisher 21 . This relation is 



m = 



3 tanh 2 (3/3 J) - 1 
2 tanh 3 (3/3 J) 



5/16 



(17) 



In table 1 we present the energy gaps of various 2 — D lattices. For sake of comparison the 3-D 
simple cubic lattice is also quoted in this table. 

Also is worth to notice that for the 2 — D lattices, the order parameter goes to zero at T c as 
(T c - T) 1 / 8 14-15,17-1^ whilc it behaves as (T c - T) 5 / 16 in the 3-D lattice 20 and as (T c - T) 1 / 2 
in the BCS theory 10 (a mean field like behavior) . In figure 4 we compare the experimental data of 
Brunei et al 13 with the exact result for the order parameter of the Ising model on the triangular 
lattice 19 . We again obtain a nice fitting of the experimental data to the theoretical curve, with the 
bonus that the zero temperature gap of the theory reproduces that experimentally observed. We 
would like also to comment that the superconducting state of pairs of spinless quasiparticles was 
considered on the basis of a model hamiltonian in MFA, by Safanov 22 . One of the results obtained 
by him is that the order parameter has a steeper rise with (1 — T/T c ) in the case of parastatistics 
than in the case of Fermi statistics. 

Finally, ultrafast dynamical optical response of YBaiCu^Oi-b was recently investigated by 
Stevens et al 23 . They found that the optical response is strongly peaked at 1.5eV, and contains two 
distinct components: one with a characteristic relaxation time of ~ 5ps, and a long-lived component 
(> 10ns) which is consistent with localized quasiparticle states at Fermi energy. For the slow 
component its differential transmittance behaves with the temperature as Ar/r oc exp (— 2Ao/fcT), 
where 2A = 3.5kT c . This corresponds to a thermal activated behavior with a fixed gap 2A . 
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We can interpret this experimental result of Stevens et al 23 , taking into account that in two- 
dimensional systems the gap only appreciably deviates from its zero temperature value very close 
to T c (see figure 1), and as we have determined in the square lattice case, this value is given by 
2A = 3.52/cT c . 

The present research has been supported partly by CNPq-Brazil. 
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Lattice 


Transitions Temperature 


Zero Temperature 




(values of kT c in J ) 


Energy Gap 


Honeycomb (z = 3) 


2.104 l 17 J 


2.85kT c 


Square (z = 4) 


2.27 I 14 l 


3.52fcT c 


Triangular (z = 6) 


3.64 I 19 l 


3.30fcT c 


Simple Cubic (z = 6) 


4.51 I 21 l 


2.66fcT c 



Table 1 : Zero temperature energy gap evaluated in this work for various 2 — 1? lattices and for 
the simple cubic lattice (3-D). 
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Figure Captions 

Figure 1: Order parameter as a function of reduced temperature T/T c . The full line is the exact 
Ising curve for the square lattice. The points are experimental results for superconducting films 
of BiiSriC aCu20% obtained by Brunei et al 13 : open circles are obtained at T c = 87K and open 
squares at T c = 78K. The dotted line represents the BCS theory and the crosses are the MFA 
results (equation (8)). 

Figure 2: Order parameter as a function of reduced temperature T/T c . The full line is the exact 
Ising curve for the square lattice. The points are experimental results for superconducting samples 
of BiiSr-iC aCu-zO-g, obtained by Staufer et al 12 : open circles are as-grown annealed at T c = 86K, 
open squares are annealed in flowing Ar at T c = 86K and the triangles are annealed in O2 at 
T c = 79K. The samples are observed at xy polarization. The dotted line represents the BCS 
theory and the crosses are the MFA results (equation (8)). 

Figure 3: Order parameter as a function of reduced temperature T/T c . The full line is the exact 
Ising curve for the square lattice. The points are experimental results for superconducting samples 
of BiiSriCaCuiO% obtained by Staufer et al 12 : open circles are as-grown annealed at T c = 86K, 
solid circles are annealed in flowing At at T c = 86K, the triangles are annealed in 2 at T c = 79K. 
The samples are observed at xy polarization. The squares are annealed in 2 a t T c = 79K but 
observed at xx polarization. The dotted line represents the BCS theory and the crosses are the 
MFA results (equation (8)). 

Figure 4: Order parameter as a function of reduced temperature T/T c . The full line is the exact 
Ising curve for the triangular lattice. The points are experimental results for superconducting films 
of BiiSriC 'aCu20s obtained by Brunei et al 13 : open circles are obtained at T c = 87 K and open 
squares at T c = 78K. The dotted line represents the BCS theory and the crosses are the MFA 
results (equation (8)). 
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